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Abstract 

We derive and discuss the constraints induced by Poincare invariance on the form 
of the heavy-quark potential up to order 1/m 2 . We present two derivations: one uses 
general arguments directly based on the Poincare algebra and the other follows from 
an explicit calculation on the expression of the potential in terms of Wilson loops. We 
confirm relations from the literature, but also clarify the origin of a long-standing false 
statement pointed out recently. 
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1 Introduction 



The effective field theory (EFT) approach has recently clarified under which circumstances 
the heavy quark- ant iquark interaction may be described by a simple potential picture. We 
call pNRQCD (potential non-relativistic QCD) the corresponding EFT. Moreover, in this 
framework the complete 1/m 2 quark-antiquark potential has been derived [|I], [2j. The final 
expression of the potential is given in terms of Wilson loops 0. Therefore, it is appropriate 
also in the situation where the quark-antiquark interaction is non-perturbative. In this case 
the Wilson loop operators may be calculated on the lattice || or in QCD vacuum models 

With respect to the classical literature on the subject (for a review see 0) the deriva- 
tion [J], 0] presents the following improvements, i) Being done in an EFT framework it 
correctly implements the ultraviolet behaviour of QCD. This is encoded into the matching 
coefficients inherited from the matching to NRQCD (on the relation between the match- 
ing coefficients of NRQCD and the heavy-quark potential see also ||). ii) Besides the spin 
and momentum- dependent potentials, it includes the 1/m and 1/m 2 momentum and spin- 
independent potentials that have not been calculated before |§ [TO]] . Hence, the obtained 
expression is complete, in) It corrects some errors present in the previous literature, which 
have propagated in several related papers until today (we refer to for a discussion on this 
point). 

In view of these new information (and in some case corrections) with respect to the 
previous knowledge on the subject, one has to ask if other results, present in the literature 
and related to the form of the heavy-quark potential, may change. In particular in fl I 



and 1 12] exact relations among the spin and momentum-dependent potentials, respectively, 



have been derived. A more recent analysis can be found in |R|. One may wonder if the 



errors found in the previous literature and/or the new terms added to the potential spoil 
the validity of those exact relations. In the present paper we will address this question and 
derive ex novo, by general considerations on the structure of the potential and by explicit 
manipulations of the expressions given in ]I|, , the relations induced by Poincare invariance 
on the form of the potentials. We shall confirm the results found in (JTTJ, and explain 
why the errors found in the literature on the derivation of the potentials do not affect the 
derivation of the exact relations among them. Since the potentials are defined up to unitary 
transformations, we shall also show that the relations found are, indeed, invariant with 
respect to them. 

The consequences of the above analysis, i.e. the obtained exact relations among the 
heavy-quark potentials, may be relevant in several situations. For instance, they allow to 
have independent checks on analytic, and lattice computations. However, in our opinion, 
the most important application is to establish and constrain the power counting of NRQCD 
under the specific situations where this EFT may be substituted by pNRQCD. This is an 
important point not only in the study of the quarkonium spectrum, for instance by means 
of lattice NRQCD, but also in the study of the quarkonium production where the validity 
of the perturbative NRQCD power counting has been recently questioned @, |14}| . 

The paper is organized as follows. In Sec. ^| we shortly describe the physical situation 
we are considering and write down the corresponding pNRQCD Lagrangian. In Sec. we 
discuss Poincare invariance from a general point of view and derive relations among the 
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different potentials without specifying their explicit form. In Sec. |] the previous relations 
are proven to be valid by direct computation on the potentials explicitly given in Jj], ||. In 
Sec. |5| we discuss the obtained results. In an appendix we clarify in some detail the origin 
of the error in the previous literature, which was found in ||, and its relation with the 
derivation of Sec. ||. 



2 pNRQCD 

In this work we consider pNRQCD in the situation where the singlet field S, describing 
the heavy-quarkonium system, is the only available ultrasoft degree of freedom. In the 
perturbative regime this situation corresponds to considering pNRQCD at the leading order 
in the multipole expansion while in the non-perturbative one it corresponds to considering 
pNRQCD without light quarks and ultrasoft gluonic excitations under the circumstances 
discussed in [|T], The pNRQCD Lagrangian reads 

£ pN RQCD = S \ id - /l s (xi,X 2 ,Pi,p 2 ) )S, (1) 



where x 3 - and pj are the position and the momentum operators of a heavy quark of mass rrij 
and spin Sj. The operator h s can be identified with the Hamiltonian of the singlet and has 
the following structure^ up to order 1/ 



2 

m 



Mx 1 ,x 2 , Pl , P2 ) = ^- + ^ + V^ (, 
lm\ Ami 



y(l,0)( r ) y%X)t T \ y(2,0) y(0,2) y{\,\) 

+ ~ + — + —T- + ZTZT^ (2) 



TOi m 2 m{ m 2 mim 2 



V {m = \ R V™(r)} + ^^L? + V^\r) + Ytf\r)U ■ S 1? (3) 



r 

,(0,2) 



V^ 2) = I {Pi V™(r)} + + K ( °' 2) (r) - Vif\r)U ■ S 2 , (4) 

= -\ { Pl • P 2 , ^(r)} - ^|^(L a • L 2 + L 2 ■ L0 + V^\r) 

MHir^i ■ S 2 - v££(r)L 2 • S : + V^^S, • S 2 + (r)S 12 (i) , (5) 

where r = xi — x 2 , Ly = r x p., and Si 2 (r) = 12r • Si r • S 2 — 4Si • S 2 . The explicit form of 
the potentials has been determined in [[l], [2[ in terms of Wilson loop operators and will be 
the subject of Sec. [§. 



3 Poincare Invariance 

Given a specific system, the generators P of space translations, the generator H of time 
translations, the generators J of rotations and the generators K of Lorentz transformations 

1 This structure already implements some of the constraints, like translational invariance, coming from 
the following discussion on the Poincare group. 
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satisfy the well-known Poincare algebra (some original references can be found in ||15|| ) 



[Pi, Pi] 


= o, 


(6) 


[P*,#] 


= o, 


(7) 


[Ji> Pj] 


= ie ijk P k , 


(8) 




= o, 


(9) 


[Jj, Jj] 




(10) 


[Pi.K,.] 


= -idijH, 


(11) 


[#,K;] 


= -iPi, 


(12) 




= ^e ^Jfc K fc , 


(13) 




= -ie y *J fc . 


(14) 



Let us consider a two-particle system of the type presented in Sec. ||. We may identify 
the generators P and J with the total momentum and the total angular momentum of the 
system respectively: 



P 
J 



Pi + P2, 



4V y 



Xi x p x + x 2 x p 2 + Si + S 2 . 



(15) 
(16) 



These definitions fulfill Eq. (y), (§) and (pTOl) . The meaning of the conditions set by Eqs. 
(0) and (|9p is trivial. The first equation constrains if to be a function of x x — x 2 , i.e. of the 
distance r of the two particles, only. The latter constrains if to be a scalar under rotations. 
Of the remaining four commutators only three are independent |T^] The constraints that 
they put on the form of the generators is by far less trivial and will be the main subject of 
this section. 

We consider our two particle system to be described by the Lagrangian (|1|) and we identify 
the generator H with 



if(xi, x 2 , pi, p 2 ) = m 1 + m 2 + h s (xi, x 2 , p x , p 2 
It has been proved in [15] that K can be written as: 

Si x pi 



(17) 



K 



2 r 

E 

i=l L 



1 + 



rrii 



tP + U, 



where {A, B} = AB + BA, and the dots indicate higher order terms in the inverse of the 
mass. These are given in closed form in the quoted literature, but do not matter here.0 U 

2 For instance, by using the Jacobi identity, it can be proven that Eq. (|l2|) follows from Eqs. (|ll]), ( |l3| ) 
and ([l4]). However, we will not make explicit use of this fact in the following. The reason is that we shall 
analyze the Poincare algebra in the framework of a 1/m expansion. In this framework, when using H and 
K up to a fixed order in 1/m, the commutation relations are not all verified at the same order in 1/m |l7f . 

3 One may wonder about the origin of the term tP in Eq. ([l8|). It is a consequence of the identification 
of Eq. ( |l7|) and the request of Lorentz invariance for the Lagrangian (|l|) (or, which is the same, for the 
Schrodinger equation). From the commutation relation (|l^ ) it then follows that the t dependence of K is 
carried by a term — tP. 
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is a function, which, due to Eq. ([13]), has to be a vector under rotations. Moreover, it is 
related to h s via Eq. fllTl). Up to order 1/m a function U that satisfies the above relations 
is 



V mi r«2 / \ mi m2 



+^ + — + — + ..., (19) 

mi m.2 

where the functions u 0) & have to fulfill the conditions 

y(°)(r) = vW(r)+vl°\r), 
V {1 '°\r) = v^ 0) (r)+vi m (r) } (20) 
V^\r) = ^°' 1) (r)+^°' 1) (r). 

The functions TJ(°\ and U^- 1 ) are arbitrary, but K obtained from them has to fulfill 

Eqs. (O)-(O). In particular, they are vectors and commute with P, i.e. they depend on 
momenta and the relative distance r only [17]]. Moreover, U^ -* commutes with V^. 



The relevant point is the following. H is now established up to order 1/m 2 , Eqs. (|T 
and (D, and K up to order 1/m, Eqs. (jT|) and (0). Then Eqs. (0), (0) and (0) have 
to be fulfilled up to order 1/m and Eq. ( |14T) up to order l/m° (consider that K starts from 
order m). The above definitions of H and K trivially satisfy Eqs. flU]), ( jl3|) and fll4|). The 



situation is different for Eq. (|T2|). In order that this relation holds at order 1/m some non- 
trivial relations between the potentials of Eq. (|2|) have to be fulfilled. Therefore, as already 
pointed out in ||18|| , the Poincare algebra may serve to constrain the form of the potentials. 



More specifically, from Eq. (|T2"D, taken up to order 1/m, after a lengthy calculation, we 
obtain relations in the spin sector of the heavy-quark potential, 

+ (22) 

and in the spin-independent sector, 

V^\r) + V^\r) - V^\r) + = 0, (23) 

-2(Vy\r) + V$ 2 \r)) + 2V$ 1 \r) - V^(r) + r^^l = 0, (24) 

\{V^\r) - V$*(r)) + ±(V™(r) - V^r)) + ^G# 0) (r) - ^(r)) = 0. (25) 

Notably the dependence on U'- 1 ' '' and U^ 0,1 ** disappears in the final relations (for 

more details see appendix [A]). If the system is invariant under mi <-> m<i transformation 
(and considering that one of the two particles is an antiparticle), we have the additional 
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conditions 



t/(2,0) 



Mr) 



1/(0,2) 
(0,2), 



V L 1 S 2 V 



mi), 
m 2 ). 



(26) 
(27) 
(28) 



In this specific situation Eqs. 
to the relation obtained in |Tl] 



) and 
and Eqs. 



reduce, up to the NRQCD matching coefficients 
(p3|)-(|2"5|) to the relations obtained in [[[^ 



In 



both cases the proof was done by explicit calculation on the potentials expressed in terms of 
Wilson loops (see Sec. |]). A proof derived from reparameterization invariance on the form 
of the scattering amplitude can be found in |13 



We note that constraints on the potentials V^'°\ V^ ' 1 ^, V r , V32 and Vs 12 could be obtained 
in this way only by including terms of order 1/m 3 in the Hamiltonian (^). These terms are 
beyond our present knowledge. 



3.1 Unitary Transformations 



The Hamiltonian and hence the potentials are defined up to unitary transformations. One 
may wonder if this ambiguity in the definition of the potentials also affects the above exact 
relations. We will show, for the example of the unitary transformation considered in [|IJ , that 
this is not the case and that the relations derived from Poincare invariance are invariant under 
unitary transformations. 

In Jl|] the following unitary transformation has been considered: 



U = exp 



W(r), 



(29) 



which transforms h s — > h! s 
are preserved by the transformation 

1 



Pi P2 

mi m 2 . 

Wh s U. We will check explicitly that the relations fl2~ID-(p5D 



h' 



h. 



2mf 



{pi,{pi,(V P W)}} 



. Up to order 1/m 2 , h' s is given by: 
1 



2mim 2 



{pi,{p l 2 ,(v;^)}} 



r{ P J 2 ,{pi(v;^)}} + 



2m^ 2mim 2 
+momentum independent terms. 



{pi,{pi(v;^)}} 



(30) 



Let us now assume a decomposition of the form V\W^ 
momentum-dependent potentials of h' s can be written as 



fif3A(r) + 5 ij B(r). Then the 



V^\r)' 


- V {1 ' 1] 

— v L2 


[r] 


+ AA(r) 


? 




_ T/(2,0) 


[r] 


+ 2A(r) 


? 


vtr\r)> 


_ T/(0,2) 


V) 


+ 2A(r) 


? 


v y\ r y 


- v 


J') 


- 4A(r) 


— AB{r 


v y\ r y 


_ T/(2,0) 

- V 


J') 


- 2A(r) 


- 2B{r 


vy\ r y 


T/ (0,2) 

- V 


[r] 


- 2A(r) 


- 2B(r 
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It then follows immediately that also the potentials obtained after the unitary transformation 
satisfy the relations (plD - (p5p . 

4 Relations obtained via Explicit Transformations of 
the Operators 



In this section we will prove the relations (pip- (|24D under the conditions (p6|)-( Pq) by explicit 
transformation of the potentials, 

^(r) = lim^MWa), (31) 



V?/\r) = ^- 2 iv lim -/ dt / ft*(t-0«yB(x 1 ,Ox(/E(x 1) *)» 

2r z T^co 1 J-T/2 J-T/2 



+|^-( V ^ (0) )> ( 32 ) 

ri ^ 1 r T / 2 r T / 2 

= ^fr-fey J_ T/ dt jf_ r/2 dt'(t-0{(ffB(x ll Ox J E(^t))> ) (33) 

V#' 0) (r) = -FF lim - / dt dt" (t -t"f ((gW(xi,f) gW(^,t))) c , (34) 

, n & — ^FF 1 /- T / 2 rT/2 

vS'°\r) = lim - / dt dt" (t - t" f (^E«(x 1 ,0^(xi,t)»c, (35) 

8 T-+oo J J-T/2 J-T/2 

7 T /-T/2 /-T/2 

Vp^W = -FF lim -/ dt"{t-t") 2 ((g^ u t")gW(^t))) c , (36) 

Z T^oo J J-T/2 J-T/2 

V£' 1} (r) = lim = / dt dt"(t-t"f((gW^t")gW^t))) c , (37) 

4 T^oo i J-T/2 J-T/2 

taken from Ref. [0] (changing, for further convenience, single in double integrals). In order 
to ease the reader we use the same notation as 0, i.e. the angular brackets (. . .) stand for 
the average value over the Yang-Mills action, Wn for the rectangular static Wilson loop of 
extension r x T (the time runs from —T/2 to T/2, the space coordinate from xi to X2): 

W a = P exp f-ig jf dz^A^z) J , dzM„ = dz°A - dz ■ A, (38) 

which we graphically represent in Fig. [l], and ((...)) = (. • • Wn) / (Wn)', P is the path-ordering 
operator. Moreover, we define the connected Wilson loop with Oi(ti), 2 (t 2 ) an d ^3(^3) 
operator insertions by: 

((0 1 (t 1 )0 2 (t 2 )0 3 (t 3 ))) c = ((0 1 (t 1 )0 2 (t 2 )0 3 (t 3 ))) - ((0 1 (t 1 )))((0 2 (t 2 )))((0 3 (t 3 ))) 

-((O 1 (t 1 )))((O 2 (t 2 )O 3 (t3)))c-((O 2 (t 2 )))((O 1 (t0O3(t3))) c -((O 3 (t3)))((O 1 (t 1 )O 2 (t 2 ))) c , 
«Oi(ti)0 2 (t 2 )» c = ((O x (t x )0 2 (t 2 ))) - ((O x (t x )))((0 2 (t 2 ))). 

The operators E i = F 0i and B i = e ijk F jk /2 (F^ = d^A u - d y A^ + ig[A^A v ]) are the 
chromoelectric and chromomagnetic field respectively. 
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Figure 1: Graphical representation of the static Wilson loop W n of extension r x T. 

The matching coefficients inherited from the matching between QCD and 

NRQCD. Their explicit form does not matter here. The only point relevant for our further 
discussion is that, from reparameterization invariance, one has ]19| 



cf = 2cf - 1. (39) 
4.1 Relations for spin-dependent potentials 

In this section we derive the relations (|21| ) -(|22| ) by explicit transformation of the potentials 



(p2\) and (|33f) . The proof is, partially, a simplified version of that one derived long ago by 
one of the authors [[11], |[]. 
We start from the identity 

«sB( Xl ,t")»=0, (40) 

which follows from parity. Then we apply an infinitesimal Lorentz boost with velocity v, 
and the point (x^t") as the origin in both systems, i.e. up to order 0(v) 

X — Xi = x' — Xi + (t' — t")\, 

(41) 

^ t - t " = t -t" + -*.[)■ v. 

The Wilson loop transforms like in Fig. |2|. The time-like paths become x' x 2 {t') = X12 — (f — 

t")v, with the boundaries t^' 2 ^' = ±T/2 — (x 12 — Xjjv. To make the last contribution vanish 
it is convenient to choose v _L r in the following of this section. The spatial lines become 
x±(s) = x 2 + s(xi — x 2 ) — v(±T/2 — t"), < s < 1. In the new system Eq. (f40|) becomes 

((gB'( Xl ,t"))Y + (([v x gE'^t")}))' = 0. (42) 

The second term comes from the transformation of the magnetic field, and the prime on the 
double brackets denotes that we have to use the transformed Wilson loop. We next drop 
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the primes for the fields and consider fl4"2|) in the old system. Subtracting Eq. (HO) from Eq. 
(0), dividing by T, and integrating over t" we get at order (D(v) 

1 / T/2 dt"(^B( Xl , *")})' " ^ d^B^O)) + ^ / T/2 dt"(([v x gE(x u t")])) = 0. 

(43) 

Notice that the second term of (42) has been just rotated back because it is already of order 



0(v). Rewriting the differences in the paths of the first and second term in Eq. (43) as 
insertions of chromoelectric and chromomagnetic fields, we get, up to order 0(v) 



f j rT/2 r T/2 

v ! - dt" / dt(t~ t" 

I T J-T/2 J-T/2 



+-e^V I ds I ' dt" 

1 JO J-T/2 



T/2 



((gB k (^t")gE\^t))) - (^B fe ( Xl , t") gE'fa, *))> 



t" + |) ((gB k (^ t") gB l {* 2 + sv, -T/2))) 
T 



t"-^)((gB l ^ 2 + sr,T/2)gB k ^ 1 ,t"))) 



T/2 
-T/2 



dt"((gW(x 1 ,t")))\ = 0. 



(44) 



I 

2 



1 

2 




gB (x„t") 



Figure 2: Difference between the boosted expression j^j ) and the original one. The arrows 
indicate the orientation of the loops, which arise after expanding up to first order of v. The 
corners of the old and the new loop have been connected by dashed lines, in order to visualize 
the form of the infinitesimal loops, which correspond to the insertions of chromoelectric 
(time-like triangles at x x and x 2y ) and chromomagnetic fields (space-like parallelograms at 
±T/2). 



Now we consider the large T limit of Eq. (44). In it has been proved that for T — ► oo 

(45) 



T J-T/2 

Moreover, we have for T — > oo 

i r T / 2 



1 r T / 2 r 1 d / 1 \ 

dt"((gW(^,t"))) = ~V^(r) + O (±) . 



r dr 



T J-T/2 



dt" 



T 



((gB k (^,t")gB l (^ 2 + sv,-T/2))) 
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((^(x 2 + S r,T/2)^B fc ( Xl ,t"))) 
v fl w b{,(°)(n| g B fc (x 1 )|0)(°) fl\ 

n^O U -fcn — £/ V J 7 

where |n)^ are the eigenstates of the static NRQCD Hamiltonian, a n is the projection on 
\n)^ of the state made by two static quarks connected by a straight string, and of the 
state made by two static quarks and the chromomagnetic field (?B'(x 2 + sr, —T/2) connected 
by a straight string (for more specification on this way of analyzing the large T limit of Wilson 
loops we refer to [jjj ||). Since, due to the different quantum numbers, this latter state has 
no overlap with the ground state \0)^°\ which we identify with the color singlet quarkonium 
state S introduced in Sec. 0, we have = 0. Therefore, the left-hand side of Eq. (f46|). 
vanishes in the large T limit. 

From the above analysis it follows that Eq. ([SID reduces, in the large T limit, to the 
identity: 



/-T/2 /-T/2 

lim - / dt" / dt(t- t") 

T^oo T J -T/2 J —T/2 



((^B(x 1 ,t ,/ ) x 0E( Xl ,t)» - ((^B(x 1 ,t ,/ ) x #E(x 2 ,t)» 



+2-V»(r) = 0. (47) 
r dr 



This corresponds to the relation^, — V{ = V^', given first in Ref. []TT]. Moreover, using 
Eq. (H), we obtain Eqs. (0)-©. 

We would like to remark that this result depends crucially on the fact that the contri- 
butions from the shifts of the space-like lines of the Wilson loop in Eq. (146]) vanish in the 
large T limit. These contributions were apparently never considered in the literature so far 
and it may be regarded as a largely lucky circumstance that they actually do not affect the 
identities involving the spin- dependent potentials. From a more technical point of view, this 
is related to the fact that the Wilson loop operators on the left-hand side of Eq. ( |46l) are 
multiplied by the factors (t" + T/2) and (t" — T/2) respectively, which are identical with 
the time differences of the correlation functions. In different situations the corresponding 
end-point contributions cannot be neglected. This fact was generally overlooked in the past 
literature and led to wrong statements and contradictions. This will be clarified in appendix 



4 The potentials V{, V 2 \ which, up to missing Wilson coefficients not considered in the early literature, 
correspond to the spin-orbit potentials, are defined as 

T j 2 T /2 

V{(r) = 11- lim i / dt f dt"(t-t")((gB( Xl ,t")xgE( Xl> t))), 

2 r T^oo 1 J_ T /2 J -T/2 
T j 2 T /2 

Y»(r) = l-. liml/ dt f dt"(t-t")((gB( Xl ,t")xgE(x 2> t))). 

2 r T^oo L J_ T / 2 J-T/2 

See also appendix |b]. 
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4.2 Relations for spin-independent potentials 



In this section we derive the relations (|23|)-(p4|) by explicit computation on the potentials 
(j34l)- (|37l) . It is convenient, for our purpouses, to introduce the following definitions^ jlO| : 

gjk = 5 iky b+ l_ _\y c 

i 1 r T / 2 r T / 2 

= lT limJ dt> «ft ff (f-O a <^(x a ,0^E*(xi,O»« (48) 

Z T->oc 1 J-T/2 J-T/2 



T jk = V k V d +r—- L -Ar\V e 



~3 r 2 



i 1 rT/2 r T/2 

= - lim - / dt' dt"(t' - t") 2 ((^( Xl) ^E fc (x 1; f)»c (49) 

The above formulae involve insertions of two chromoelectric fields, together with quadratic 
factors in the time differences. Therefore, we have to consider the Lorentz transformation 
up to second order in the velocity v. 

We start with an expression like fl4"D|), but with a chromoelectric field instead of the 
chromomagnetic one, i.e. 

((i<7V.E( Xl ,0)>. (50) 

The velocity v will be used in the Lorentz boost. The longitudinal component of the chro- 
moelectric field which appears in (|50|) is not changed under this transformation, this greatly 
simplifies the derivation. Another reason for taking the scalar product with v will become 
obvious below. The Lorentz transformation with (xi, t") as origin is now needed up to order 
£>(v 2 ). It reads 



X - Xi = x' - Xl + (t' - *")V + -[(X' - Xl ) • V]V, 



v 2 ' 



(51) 



t-t"= + (f - 1") + (x - Xl ) • v. 
In the new system the time-like lines of the Wilson loop become 



with the boundaries 



< 2 (0 = - («' - Ov - j[(%,2 - x,) ■ v]v, (52) 



i^'-ifl + yU-lxu-xO-v-fi". (53) 



In terms of the potentials defined at the beginning of Sec. H, we have 



\, « -\v„ v^ = v e , v^ = -v h + z -v C) v^ = -v c . 
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The spatial lines become 



1 



x' ± (s) = x 2 + s( Xl - x 2 ) - v(±T/2 - t") + | (1 - s) [( Xl - x 2 ) • v]v, < s < 1. (54) 

We now write (|50|) in the new coordinate system. Recalling that the longitudinal component 
of E is not transformed, we obtain 

((^•E'( Xll f)}}'. (55) 

Again the prime denotes the transformed Wilson loop. We next drop the primes for the 
field in the boosted system and consider the difference between the expressions ( |5"5] ) and 
([501), divide by T, and integrate over t" from — T/2 to T/2. This has to vanish by Lorentz 
invariance: 

- fj^ dt"{(igv • E( Xl , t")))' - i df{{igv ■ E(x lf 0» = 0- (56) 

The figure looks like Fig. ^|with B( Xl ,t") replaced by E( Xl ,t"), but the shift of the Wilson 
loop lines has now to be considered up to order v 2 . 



/0' 



x + Ax 



Figure 3: Distortion of the path showing the appearance of the field insertions and the 
covariant derivatives. 

In the Abelian case the evaluation of the left-hand side of Eq. (|56| ) could be easily done 
according to the self-explaining scheme 




11 



where C is the original Wilson loop and C the boosted one. More specifically, considering 
the time-like paths only, we get J C i_ c = j E da., where the integral is now over the area 
enclosed by the time-like lines of the paths C and C . The two paths differ by order v. In 
order to get the correct result up to order v 2 one has to expand E(xi j2 — s(t' — i")v — s[v ■ 
(xi >2 - xi)]v/2,t') = E(x 1)2 ,t') - s(t' - t")(v ■ V)E(xi j2 ,f) + 0(v 2 ). Here s parameterizes 
the transversal direction between the lines C and C; it has to be integrated from to 1, 
which gives a factor 1/2 in the second term. Analogously we may calculate the contributions 
coming from the shifts in the spatial lines (these, however, involve also chromomagnetic 
fields). 

In the non-Abelian case one can apply a method similar to that used in the derivation of 
the non-Abelian Stokes theorem ||20|| . Distort the line at, say, x + Ax by adding appendices 



as shown in Fig. |3|. This gives a path at x, together with insertions of chromoelectric fields, 
represented by the paths around the small rectangles. They start and end in the middle of 
the rectangles and give the chromoelectric field at the position x + Ax/2. Together with the 
small lines connecting the midpoint with the line at x one has 

e - l ,A(x).Ax/2 E fe ^ + ^ e , 9 A(x).Ax/2 = E * (x) + Ax {(y^k^ _ 

^[A(x),E fc (x)])+... = E fc (x) + ^-[D x ,E fc (x)] + ..., (58) 

with the covariant derivative D x = V x — igA(x). 

Applying the above considerations to the left-hand side of Eq. (56) one finds that the 
contributions of order v 2 vanish by time reversal invariance, because the integrals contain 
factors (f — t"), which are linear in the time variables. The spatial shifts do not contribute 
in the limit T —>■ oo. This can be seen explicitly by using the same methods and arguments 
of Sec. [£T[ The remaining contributions of order v 3 are, in the large T limit, 

1 r T/2 r T/2 

= ± lim -i- / dt> / dt"{t'-t"f 



T^oo 2T J-T/2 J-T/2 

x (([v ■ D X1 2 , ■ E(x li2 , if)] gv ■ E(x x , t"))) c , (59) 
E® = - lim — (r • v) / & / dt"((gv ■ E(x 2 , t') gv ■ E(x 1; t"))) c , (60) 

T^oo 21 J-T/2 J-T/2 

EF> = E® + E®-2E&\ (61) 
with 

,^ i r T / 2 r T / 2 r T / 2 

E\J = - lim — / dill dt" df'if -t"){t'" -t") 

Ok T^oo2T J-T/2 J-T/2 J-T/2 V A ' 

x ((gv ■ E(x„ If) gv ■ E(x fc , if") gv ■ E( Xl) t"))) c . (62) 

The terms E^- 1,2 ^ arise from the first order shift (with respect to v) of the paths, together 
with the first order expansion of the chromoelectric field. E^ belongs to the second order 
shift in ([52D at the line x 2 , while E^ is due to the square of the integral in (|57D. 

We next rewrite the commutator of the covariant derivative in E^- 1 ' 2 '. This is rather simple 
for E^ 2 ': the loop with the insertion of [D£ , E fc (x 2 , t')\ may be written as V X2 acting on the 
loop with the insertion of E fc (x 2 , t'), minus V X2 acting on the loop alone. The last derivative 
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can, in turn, be expressed, in the large T limit, by an additional insertion of a chromoelectric 
field igE- J (x 2 , t'") and integration over t'". In this way becomes a sum of a derivative of 
a correlator with two fields and of a correlator with three fields, = + , with 

= - lim -(v-V X2 ) / dt' dt"{t'-t"f 

T^oo 21 J-T/2 J-T/2 

x((( 7 vE(x 2 ,0^v-E(x 1 ,t"))) c , (63) 

/rt\ j rT/2 rT/2 r T/2 

EP = lim — / dt' / dt" / dt"'(t' - t"f 

6 T^oo2T J-T/2 J-T/2 J-T/2 V ' 

x ((gv ■ E(x 2 , t') gv ■ E(x 2 , t'") gv • E(x ls *")»c (64) 

For E^ 1 ) the situation is slightly different, because V£ acting on the whole loop together 
with the insertions also differentiates E(xi, t"). The fields v • E(xi, t') and v • E(xi, t") enter 
symmetrically, therefore differentiation gives twice the same expression. Instead of (|63|) we 
thus obtain E<& = E^ + with 

1 r T / 2 r T / 2 

Eg* = lim— (vV Xl )/ dt' rft"(t'-t") 2 ((^v-E(x 1 ,0^v-E(x 1 ,t"))) c ,(65) 

T^oo 4i J-T/2 J-T/2 

r T/2 r T/2 rT/2 

- / dt' dt" / dt"'(t' - t"f 

T-^oo 4T J-T/2 J-T/2 J-T/2 

x ((gv • E(x 1; t') gv ■ E( Xl , t"') gv ■ E(x ls t")»c (66) 
The correlators i?! and can be directly written in terms of the expressions S^ k and T jfc 



= lim 



in (|48|)-(|49|). As intermediate result we obtain 

Ef ] + Ei 1] = i(v ■ V r )\*(S ik - T jk )v k . (67) 

The contribution E^ can be expressed by the static potential if we rewrite the insertions 
of the chromoelectric fields as derivatives acting on the loop. It can be written in the form 

£(3) = -l( r ■ v )(v • V r )^— ^V (0)/ (r) = ~(v • V r )v J ' ( 5 jk V^(r) - —V^'(r)\ v k . (68) 
2 r 2 [ r J 

We next consider the expressions which contain insertions with three fields, namely E^ , E^ 

n(4) 



and the three terms in E^ . In i% we can symmetrize the factor in front with respect to t! 
and t", i.e. replace (tf - t")(f" - t") -> (tf - t"f/2. Then one has + E$ = 0. In E® 
and ^22^ we exchange the coordinates xi and x 2 by applying CP, furthermore we exchange 
the integration variables t" and t'" . The remaining three terms then all involve the same 
correlator. Symmetrising with respect to t',t" we get 

4 2) + E® - 2Ej 4) 



22 ^-^12 

3z 



[T/2 


fT/2 


/•r/2 


/ dt' 


/ dt" 




l-T/2 


'-T/2 


J-T/2 


x {{gv ■ 


■E( Xl ,t') 


\gv ■ 



lim 



3 , „ , r T / 2 , , r T / 2 

hm — (v • V 



lI/\2 



i-rri i-i2 

r - ..... // W ^(f-O^'ECx^Oyv.ECxi.O^ 

T^oo 4J J-T/2 J-T/2 

3i(v • V r )v J T J V\ (69) 
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In the third line we wrote the insertion J ■ ■ -igv ■ E(x2,t") dt" as differentiation of the loop 
with respect to x 2 , in the last line we introduced the definition of T jk in fl49|) . Lorentz 
invariance requires that the sum of (p7|), (pE|), QS5I ) vanishes. This gives 

S jk + 2T jk + t_ y(o) _ LiLy(°) '( r ) v fc = 0. (70) 
2 2r j 

The term in curly brackets is just the expression which has to vanish according to Eqs. 
(HD-flU) and 0. 

Eq. ([T(]) is only slightly weaker than the statement of the vanishing of the curly bracket. 
Let us introduce the combinations 

f(r) = V d (r) + ~V b (r) + \v^\r) - ^V (0)/ (r), g(r) = V e (r) + \v c {r) + ~l/ (0) '(r). (71) 



These are the functions that have to vanish according to Eqs. (p3|)- (p4|) and [12| if S 3 and 
T jk are decomposed into V&(r), • • • , V e {r). In our case (^) gives two differential equations 
which appear as factors of 4iv 2 (r • v)/r 2 and 4i(r • v) 3 /r 4 , respectively. They read 

rf (r)/2-g(r) +rg'(r)/Q = 0, g(r) - rg'{r)/2 = 0. (72) 

The general solution is g(r) = ar 2 , f(r) = 2ar 2 /3 + b with a,b = constant. Only a = 
appears physically reasonable, while the constant b may be reabsorbed into the static 
potential. 



5 Conclusions 

In this work we have derived identities among the potentials of pNRQCD (as given in Sec. 
H) either by direct application of the Poincare algebra to the potentials in their general 
form or by explicit Lorentz transformation of the potentials in terms of Wilson loops. In 
this way we have proved that the identities found long ago for the spin-dependent |T 



and spin- in dependent potentials |T2| are correct, despite several subtleties involved in the 
manipulation of Wilson loops overlooked by the literature (as first pointed out in and 
discussed in appendix ||). We have thereby also proved that Lorentz invariance is the reason 
behind these two set of identities, providing a unified framework. 

From the general arguments of Sec. |3| the identities (|2~ID -(|25D are the only ones that can 
be derived, between the potentials up to order 1/m 2 . The identities that could constrain the 
potentials of order 1/m and the momentum independent 1/m 2 potentials would involve 1/m 3 
potentials, which are, at present, unknown. In this respect the "brute force" method used in 
Sec. |] may be more useful. It may be used in order to constrain Wilson loop operators that 
are relevant for the calculation of the potentials up to order 1/m 2 , but do not necessarily 
show up as potentials themselves. In particular, one could apply transformations similar to 



those discussed in Sees. [O] and fO| on Wilson loops with two or more field insertions. This 
would result in constraints involving Wilson loops with three or more field insertions. These 
type of Wilson loops are known to contribute to the momentum-independent potentials of 
order 1/m 2 0. 
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The complementarity of the two methods discussed in this work, which, however, are 
both applications of Poincare invariance, can be made clear by noticing that the identities 
(PT|) and ([22D proved in Sec. [3] actually differ from the identity ( pETp proved in Sec. ( |4.1| ). 



Indeed, combining Eq. P7| ) with Eqs. (PT|)-(p2|) we get Cg = 2c^p — 1, which, therefore, may 
be derived in this way from Poincare symmetry without using reparameterization invariance 
arguments. 

Finally, we comment on the power-counting issue mentioned in the introduction. Eqs. 
(PT|)-(P^) contain combinations of 1/m 2 potentials and the static potential. Since scales 
like mv 2 in the heavy quark velocity v, as usual in non-relativistic bound states, then the 
combinations of potentials of dimension three appearing in Eqs. (pll)- (|22|) must scale as 
m 3 v 4 and those of dimension one appearing in Eqs. (p3|)-(|24]) as mv 2 . This result is not 
trivial. It implies that the considered potentials are suppressed by an extra power of v with 
respect to the "natural" power counting based on their dimension. In perturbative QCD this 
extra suppression factor is typically provided by the coupling constant, a s ~ v. Poincare 
invariance tells that an extra suppression factor v must be dynamically generated also in the 
non-perturbative regime. 

Acknowledgements. N.B. and A.V. thank the Institut fur Theoretische Physik of 
the University of Heidelberg for the warm hospitality and the Alexander von Humboldt 
Foundation for support during the first stage of this work. 
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A Potentials and Lorentz Transformations 

In this section we show in more detail how the dependence on Uw, U'- 1 ' -' and U^ ' 1 ^ cancels 
in the final relations (PT|)-(25). In order to be definite we will assume 

U (0) (r) = U W (r)r, (73) 
U (1 ' 0) (r, P i) = {s(r)^Pl} + {/(r), Pl } + Mr)r, (74) 

U^(r,p 2 ) = {g( r )-,pl) + {f(r),p 2 } + h(r)v, (75) 



where g, f and h are arbitrary functions. From Eq. ( II ) it follows that one must have 
g(r) = 0. Inserting Eqs. (|I"9|), (pT8|), (|P7D in Eq. (|12|) we get, up to order 1/m, the following 
five equations: 

r^vf\r) + 2Vl 2 2 ' 0) (r) - V^/Ur) - r^-U {0) (r) = 0, (76) 
dr dr 

m<°> (r) + 2 V$ 2) (r) - V$ X) (r) + r^-U {0) (r) = 0, (77) 

-v^ (r) - 2^' 0) (r ) + v£ A) (r) - 2 V$ 0) (r) + (r ) - U (0) (r) = 0, (78) 
-^ 0) (r) - 2 V^' 2) (r ) + V$ 1] (r) - 2y^°' 2) (r) + (r ) + (r) = 0, (79) 

£Wly(o)( r ) = . (80) 
r ar 

Adding Eq. (|7§) to Eq. (|77|) and Eq. ([78]) to Eq. (|7p]) we get Eqs. (|23|) and (£§. Subtracting 
them and combining the result, we get Eq. fl25|). Eq. (|^) gives 

f(r) = 0, (81) 

therefore it follows that U*- 1,0 - 1 and U^ 0,1 ^ do not depend linearly on pi and p 2 respectively. 
The above derivation may be extended to more general functions U^ 1 ' - 1 and U^ ' 1 - 1 that admit 
a power expansion in the momenta. 



B Subtleties and Incorrect Proofs 

In this section we clarify a false statement, due to a wrong treatment of end-point string 
contributions, which can be found in the literature. Let us consider the function |§ |nj |Tl 



7 r k 1 r T / 2 r T / 2 

Viir) = -e w '- lim - / dt dt\t-t")(gW{x u t")gW{^t))) ) (82) 



related to the spin-orbit potentials introduced in Sec. |] through the relation c^V^ = fV^'^. 
In the review || one finds the statement that the factor t" in the integrand may be dropped, 
i.e. that the value of V{ does not change if we replace (t — t") with t. Indeed, in this second 
form the potential V{ can be found, for instance, in ||, p70|| . It has been pointed out and 
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Figure 4: Difference between the expression fcjQ ) shifted accordingly to x 2 — > x 2 + Ax 2 and 
the original one. The arrows indicate the orientation of the loops, which arise after expanding 
up to first order o/Ax 2 . The corners of the old and the new loop have been connected by 
dashed lines, in order to visualize the form of the infinitesimal loops, which correspond to 
the insertions of chromoelectric (time-like rectangle) and chromomagnetic fields (space-like 
triangles at ±T/2). 



proved in M that this statement is wrong even at leading order in perturbation theory. In 
the following we will make clear how it was possible to arrive at that erroneous statement. 

Let us start from the identity (0). Applying the differential operator V^ 2 to it and 
writing the differentiation of the time-like lines as insertions of chromoelectric fields and the 
differentiation of the space-like lines as insertions of chromomagnetic fields, we obtain 

T/2 

dt((gB l (x u t")gW(x 2 ,t))) 

-T/2 

-ie jkl r k f 1 ds (1 - s) ((#B'(x 2 + sr, T/2) gBHx^t"))) 
Jo 

+ie jkl r k ^ ds (1 - s) ((gB\^,t") gB l {x 2 + sr, -T/2))) = 0. (83) 
Jo 

Assuming that the last two lines, once inserted in Eq. (^3), vanish in the large T limit (as 
tacitly assumed by previous authors), leads immediately to the conclusion that the factor t" 
may, indeed, be dropped in the integrand of flS2|). 

The key point is that, as argued in 0, the end-point string contributions, corresponding 
to the small triangles at time ±T/2 of Fig. |], do not vanish in the large T limit in Eq. (|82|) . 
This is easy to check at lowest order of perturbation theory, where we have: 



: / dt ((gBUxxX) sE J (x 2 , *))) P =, - 2iX^e i3l v l 

J -T/2 7T 

/ 1 1 \ 

,(84) 



\{t" - T/2) 2 - r 2 - ze] 2 \{t" + T/2) 2 - r 2 - te 
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ie jkl r k C ds (1 - s) ((gB*(xi, t") gB l (x 2 + sr, -T/2)» P = t - 2i^^e ijl r l 

JO 7T 

X [(t" + T/2) 2 -r 2 -ze] 2 ' ^ 
ie jkl r k f 1 ds (1 - s) ({gB l (x 2 + sr, T/2) gBH^, if'))) P = t - 2i < ^e ijl r l 

JO 7T 

1 

\{t" - T/2) 2 - r 2 - ie] 



x iu„ 7F7^ ^ ,, 12 » ( 86 ) 



with C/ = (iV 2 - 1)/2JV C = 4/3 in QCD. It is clear that all three lines of Eq. (|3|) give 
the same kind of contribution in leading order perturbation theory and that the identity is 
satisfied only if all three pieces are taken into account. 

The argument can be extended to any order of perturbation theory. If we analyze the 
above end-point string contributions by decomposing them into the eigenstates \n)^ of the 
static NRQCD Hamiltonian, we obtain for T — > oo 



i 



/■T/2 r 

/ dt"t" ({gB^ h f)gB\^ + sv,-T/2)))-{{gB\^ 2 + sv,T/2)gB\ Xl ,f))) 

1 J -T/2 



1 a^W^l^xQIO)^)-^^^)^!^^^)^)^) /IX 

^) — \t r ( } 

U n ^o tin — -C/ V 1 7 



with the same notation as in Sec. |4.1| . The ground state lO)^ -* has no overlap with the state 



made from two static quarks and a chromomagnetic field gB*(x 2 + sr, —T/2) connected by 
a straight string, hence = 0. However, contrary to Eq. (|46D , there are still contributions 
that do not vanish for T —>■ oo. This general argument definitely proves that the end-point 
string contributions of Eq. (p3|) do not vanish and, therefore, we cannot drop the factor t" 
in Eq. 

So, we finally find that a forgotten contribution corrects a wrong proof where necessary, 
but cancels in a well-known relation (Eq. (pTTD), where it would be very disturbing. Future 
manipulations of the above type should, therefore, be done with the appropriate care, and 
old ones revisited appropriately. 
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